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ABSTRACT 



This paper presents an independent verification of the 
collocation method as a technique for calculating the lift 
on an oscillating airfoil in an unstaggered cascade 
immersed in transonic flow. This method was originally 
proposed by Gorelov. Results presented here differ some- 
what from those presented by him. Two formulations are 
shown; one is purely numerical, the second employs an 
analytic expansion for small frequency. 
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LIST OF SYMBOLS 



local speed of sound 


II 


speed of sound in the uniform flow 
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blade semichord 


III 


elementary function used in 
collocation solution 


V 


specific energy ("head") 


II 


function describing the surface of the 
airfoil as a function of time 


III 


function specifying location of blade 
surface in the vertical axis 
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Mach number = -* — ‘- 
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II 
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x-component of velocity 


II, III 
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respectively, solved so as to satisfy 
the tangential flow conditions 
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Ill 



small disturbance velocity 


III 


y-component of velocity 


II, III, IV, V 


general velocity vector 


II 


small disturbance velocity 


III 


vertical velocities due to the reference 
and adjacent blades respectively, 
determined from the tangential flow 
condition 


V 


(p^/ a constant used in Gorelov's 
approximation of the transonic flow 
potential 

horizontal coordinate, may be 
non-dimensional 


IV, V 


mp (transformed interblade distance 
in Gorelov's approximation) 


V 


blade leading edge 


IV 


center of pitch of the unstaggered 
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vertical coordinate, may be 
non-dimensional 


IV,V 


vertical coordinates attached to the 
reference and adjacent blades 
respectively, may be non-dimensional 


IV, V 


transformed vertical coordinates used 
in Gorelov's approximation, attached 
to the reference and adjacent blades 
respectively. z = my, = my^ 


IV, V 


substantial derivative w.r.t. time 

_ 3 , ax 3 , ay a 

at at 3x at ay 

= - — + U - — + V — 

at 3 x ay 


II, III 


2 

"of the order of magnitude of o) " 


V 


angle of attack 


II, III, IV, V 


maximum amplitude of pitch- oscillations 


IV 
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V 



ratio of specific heats, c /c 

p V 

Dirac Delta function 

_ 1 when i = 0 
0 when i 0 

cos (-x) 


V 

V 


-1/, 

cos (1-X*) 


V 
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cos (-s) 
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V 


interference potential coefficients 
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V 


k/ra^ 


IV, V 


Fourier coefficients describing the 
motion of the reference and 
adjacent blades respectively 


V 


angular frequency of oscillation 


IV, V 


density 


II 


phase angle 


V 


. . 2 
cascade solidity, — 


V 


general velocity potential 
uniform flow velocity potential 


II, III, VI 


steady flow perturbation potential 


III, IV 


perturbation potential in collocation 
solution due to reference and 
adjacent blades respectively 


H 

> 

> 


tranformed potentials 


V 


oscillatory flow potential 


III, IV, V 


transformed oscillatory potential in 
Gorelov's coordinates corresponding to ip 


V 


interference potentials due to 
reference and adjacent blades 
respectively 


V,VI 
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V 



transformed potentials in Gorelov's 
coordinates, coresponding to and 



k(l-m) 



m 



w + 3 



Gradient operator. 



1 and 3 



3x ‘ 9y 

are unit vectors in the x and y direction 
respectively 
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Computer Variables 



DK 

DLAMDA 

DM2 

DR 

ETA 

ETASTR 

IPT 

N 

NF 

OFFSET 

OMEGA 

QALPHA 

QCONST 

QDCL 

QDCM 

QDK 

QEXP 

QINTAP,QINTRP = 

Q1ABCF,Q1RBCF = 
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Q1RBP,Q1ABP 
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Reduced Frequency, k 

X 

2 

m 

mp = X* 
n 

Print Parameter 
n 

not used 
r 
0) 

0 + i(X-k) 
e^^ 

Cl 

a 

Cm 

a 

0 + ik 
0 + ix 

Variables used to transmit boundary 
condition integrals 

Interference coefficients for adjacent 
and reference blades 

Right hand side vector in collocation solution 

Known matrix of integrals in collocation 
solution 

Not used 

Not used 

^-ixx 

e 



11 



Q2PT 

RHO 

RHP 

SIGMA 

TAU 

XASTN 

XSTN 



= Not used 

= vertical displacement 
= local input variable for rho 
= a 

= T 

= Current x station in adjacent blade 
coordinates 

= Current x station in reference blade 
coordinates 
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I . INTRODUCTION 



The analysis of unsteady transonic flows in aircraft 
turbopropulsion is an area of intense current interest. 

Rising fuel prices and increasing thrust requirements both 
point toward the need of turbomachinery capable of performing 
well with transonic or supersonic internal flow. But, 
increased flow has increased both the costs and uncertainties 
of engine designs. Flutter problems have already become a 
major consideration in engine development. Problems unfore- 
seen in earlier days of turbine engine production have 
caused long development delays, or forced acceptance of 
engines producing less than their initial design thrust. 

These uncertainties cannot be avoided when an attempt is 
made to extend the state of the art, but they can be reduced 
by extending the range of analytical modeling. 

Such extension must now be done piecemeal. The three- 
dimensional flows in turbomachinery, including the simul- 
taneous effects of boundary layers, rotation, finite blade 
thickness, spanwise Mach distributions, and shocks, are 
well beyond present capability. Perhaps one day complete 
analysis will be practical, but it is not today. The best 
that can be done now is to approach the problem from one 
aspect at a time. Flow through a two dimensional cascade 
has been a useful tool in this partial analysis. 
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This thesis was originally to have been an extension of 
the work of Elder [1] and Schlein [2] to the case of a 
staggered cascade. Their work, based on Teipel's [3] 
linearization of the unsteady transonic small perturbation 
equation, analyzed transonic flow through oscillating 
unstaggered cascades by use of the collocation method. 

While the problem was easy to state, it was difficult to 
solve. Both Elder and Schlein had encountered difficulty 
in employing the collocation method. Therefore, it was 
decided that verification of the basic collocation solution 
presented by Gorelov [4] using a different linearization 
would-be a worthwhile goal in itself. 

The following investigation presents a verification 
of the development in [4] , along with numerical results and 
suggestions for further work. 
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II. UNSTEADY TRANSONIC FLOW THEORY 



Considering inviscid flow only, the following four 
equations govern the aerodynamic flow problem at hand; 
The equation of state 



p = pRT 



(II-l) 



and the equations for the conservation of 



1. Mass : div(pv) + 



i£ 

9t 



= 0 



(II-2) 



2 . 



Momentum: 



Dv 1 
Dt p 



Vp 



0 



(II-3) 



3. 



t:. DS 

Energy: — 



0 



(II-4) 



where 

V = velocity 
p = pressure 
S = entropy 

R = universal gas constant 
T = temperature 
t = time 

p = density 

D _ ^ ^ ^ ^ ^ 

Dt ” at 9x 3t 3y at 

1 >h 0 analysis starts with 'a uniform flow from infinity. 
This flow has velocity parallel to the x-axis. This 
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1 




formulation can be simplified by working with the total 
velocity potential, where 

3 $ 

u = ~ ~ ^ component of velocity = (II-5) 

V = = y component of velocity = (II-6) 

Thus, the initial uniform flow is represented by the uniform 
flow potential 

<J> = U X (II-7) 

o o 



This notation may be applied to the conservation equations 
for mass and momentum. The equation for conservation of 
mass 



div(pv) 



(II-8) 



becomes for two-dimensional unsteady flow 



9 (pu) 3 (pv) I 3p _ 

3x 3y 3t 



+ U |§ + V |£] + P(|^ + 



3t 



3x 



3y 



— ] 

3yJ 



0 



(II-9) 



but 



i£ 

3t 



+ 



£ 

3x 



V 



l£ 

3y 



= ^ 



Dt 
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and 



Thus 



and 



u = $ and V = $ 



^ + p ( + <5 ) = 0 

Dt XX yy 



$ + (|i 

XX yy 



1 Dp 
“ p Dt 



(II-IO) 



The speed of sound is given by 



2 _ dp 

' dp 



Thus 



Dp _ ^ ^ Dp _ Dp 

Dt dp ' Dt 2 Dt 



Applying this to equation (II-IO) yields 



$ + $ 

XX yy 



1 Dp 
2 Dt 



pa 



(Il-lla) 



$ + $ 

XX yy 



pa 



•(UP^ + vp + p^) 



(Il-llb) 



pa 



--jl (V$) • (Vp) + p^] (II-llc) 
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p = 
X 



9x 



p = i£ F = ^ 

y 3y ' t 3t 



where V is the gradient operator. 

Laying this aside for the moment, consider the momentum 
equation (II- 3) 



Thus 



Dv = 3v rt\ 

Di ^ + (V • V) V 



and 



+ (V 



V) v] = Vp 



( 11 - 12 ) 



so 



= V$ 






= V 



_ 3 £ 

3t 



(11-13) 



(V • V) V = 



V ^ - V X ( V X v) 



where 



2 2 2 
V = u + V 



= V * V 
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For irrotational flow 



V X V = 0 



thus 



(v • V) V = (11-14) 



Thus 

2 

^ + V[<J)^ + ^] = 0 (11-15) 

which after integration along a streamline becomes 



/^ + <l>t + ^ = F(t) (11-16) 



For uniform flow from infinity F(t) 
the final result is 




and thus 



p 



+ 4 >. 





(11-17) 



Differentiation with respect to t gives 



20 . 



(11-18) 



- p($ + 1 iZl) 

^ ^ tt 2 at ^ 



From (II-3) 



Dv 

■vp = p 5^ 



(11-19) 



Substitute (11-18) and (11-19) into (II-llc) to obtain 



p + <i> 

XX yy 



1 r . Dv . ^1 av^. 

^[(V») 5t «>tt + I 

ci 



( 11 - 20 ) 



This may be further simplified 



Dv 

Dt 



av 

at 



+ (v • V) V 



Hence 






$ + <5 

XX yy 



(V<f>) • (|^ + jVV^; 



^ *tt * I H-> 



Expanding terms 



(V») • (|f) 



= V<|) • (^V$) = + 

9t X xt y yt 



v<j> 



vv 



2 2 

(|)(|) 

X XX y yy + x y xy y X yx 
2 2 2 2 






$ 



2^ ^ ^ 

X XX y yy + ^ y xy 

2 2 2 
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1 1 

2 ^ ' 2 • <''»>) 



— t" [ $ $ . + $ , $ +$ $ . + <J> $ ] 

2 X xt xt X y yt yt y* 



*x*xt ■" *y*yt 



The final result obtained by combining terms is 



$ $ 2$ $ $ 
(1 - - 4 ^) + (1 - -^) ^ y 

^ 2 ^ XX 2 ’ yY 2 

a a a 



2$ 



2 $ 



^ <l> . - i , = 0 



2 xt 2 yt 2 tt 

3 . 3 3 



( 11 - 22 ) 



This result is valid for irrotational , inviscid, two- 
dimensional, unsteady, compressible flows where gravity has 
been neglected. 
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III. SMALL PERTURBATION THEORY OF TRANSONIC FLOW 



A. GENERAL CASE 

A thin body at a small angle of attack will cause only 
a slight disturbance in the fluid. A flat plate is an 
example. Consider flow past a flat plate at angle of 
attack', a . 




The flow at (2) must be parallel to the plate. To achieve 
this, small disturbance velocities u' and v' must be added 
to the free stream velocity yielding 

u = U + u' 
o 

V = V ' 

The potential of the disturbed flow may be considered 
as the sum of the uniform flow potential, = U^x, and the 
disturbance potential, <J> 



$ 



$ + (f) 

o ^ 



(III-l) 
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Thus 






X 



U + u' 
o 



= V 



If <j> is a function of time,, then 



(III-2a) 



(III-2b) 



\ 



This result may be substituted into (11-22) leading to 



(U + u')^ 2 (U + u')v 

[1 - ^ + [1 - - 2 ♦ 

a a a 



xy 



(U^ + u’) 

2-2-= 4, 



2v 



xt ,2 '♦'yt ,2 ^tt 

3 > 3 . 



4)^^ = 0 (III-3) 



This expands to yield 



(U^v + u'v ) 



U ^ + 2U u' + u'^ „2 .. 

[1 2 ] 4 , + [1 _ 2 _] 4 ) - 2 — 2 ^ 4 ) 

2 ^xx ^ YY 2 ^xy 

3 3 3 



2v 



(U + u') 

2 « (p — (p - ~ 0 

2 ^xt 2 ^yt 2 ^tt 

a a a 



(III-4) 



Equation (III-4) may be further simplified as shown by 
Landahl [ 3 ] . All non-linear terms except the 4>v'J*w Pi^oduct 
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term can be neglected yielding the following transonic 
small disturbance equation 



[(M^-1) + (y+1)M^ ^](|) 

U ^xx 
o 



— (p + 

yy 



2M , 

— 2 . 1 ^ 1 . 

% 



tt 



= 0 



(III-5) 



where a^ is the velocity of sound in the free-stream, and 
Y is the ratio of specific heats, and M’<=\Mach number. 



B. BOUNDARY CONDITION 

The tangential flow condition requires that the flow be 
tangent to the airfoil surface at each instant of time. This 
means that no fluid may flow through the surface of the air- 
foil and is expressed by the condition 



DG 

Dt 



0 on G(x,y,t) 



(III-6) 



where 



G(x,y,t) describes the surface of the body 
as a function of time. 

For a thin airfoil restricted to small oscillations, this 
may be written as 



G- = y - H(x,t) 



(III-7) 
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where 



H(x,t) is the function describing the 
position of the airfoil. 

H(x,t) can be written for harmonic pitch oscillations as 

H(x,t) = R.P. [a (x-x ) e^^^] (III-8) 

o o ' 

where the time-varying angle of attack a(t) is given by 

a(t) = R.P. 

and OL^ = maximum amplitude of pitch oscrllation 

x^ is the pitch axis 

V is the angular frequency of oscillation 
i = /-I 

R.P. = "real part of" 

Inserting (III-8) into the flow tangency condition 
(III-6) gives, after linearization, 

(j)y(x,0) = v(x,0) = + iv(x - x^) ] e^^^ (III-9) 

on y = 0 

This is a condition for the normal velocity to be prescribed 
at the airfoil's mean position y = 0. 

C. NONDIMENSIONALIZATION 

The terms in equations (III-5) and (III-9) are dimensional. 
For the following calculations it is convenient to use non- 
dimensional quantities. Define non-dimensional time and 
length to be 
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where 



X = 



X 

c 




(III-IO) 



= uniform velocity from infinity 
c = reference length (blade semichord) . 



The velocity potential in equation (III-5) may be non- 
dimensionalized as follows. Let 









Hence : 



(j) = U^C(t> 

o ^ 

A = U C^-(-) 

o c 

= u (III-ll) 

o^x 

and similarly for the other derivatives in (III— 5), yielding 
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0 



[(M^-1) + (Y+l)M?(j)-] 4)— -- 4)--+ 2M^(1)--+ 



(III-12) 



This equation is non-dimensional. 

The boundary condition given in equation (III-9) may 
be non-dimensionalized in a similar fashion 



(x, 0) 



a lu + iv (x - 
o o 






ivt 



(III-9) 



Thus 



. — c 

U C(b = a [U + ik— (cx - cx^] e ^o 
o ^y o o c o 



(III-13) 



where 



k 



^ = Strouhal number or reduced 

frequency 



U c<}>-'- 
o ^y c 



a U [1 + ik(x - 
o o 






ikt 



(III-14) 



Thus 



(P- = v(x,0) 






+ ik(x - 






ikt 



(III-15) 



Because the final operations are linear in a^, set - 1, 
yielding 

— — 1 Ic "f" 

i- = v(x,0) = [1 + ik(x - x^)] e^^^ (III-16) 
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The overbars denoting nondimensional quantities will be 
dropped from the remainder of the paper. All further 
quantities shall be assumed appropriately non-dimensional. 

This yields the following final equations 

t(M^-l) + (Y+l)M^(p^]4,^^ - = 0 

(III-17) 



and 



4>y(x,0) = v(x,0) 



1 Tc 'h' 

[1 + ik(x - Xq)] (III-18) 



where 



all quantities are nondimensional and 

a = 1 

o 

D . HARMONIC OSCILLATIONS 

In the case of harmonic oscillations, equation (III-17) 
may be simplified still further. 

Let 

-i Ir 4- 

p = !}) + R.P. [ij;e ] 

where 

(p = non-dimensional steady flow potential 
'P = non-dimensional oscillatory flow potential 
R.P. = "real part of" 
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Equation (III-17) then becomes 







X 



0 



(III-19) 



For M close to 1, this is a nonlinear mixed elliptic- 
hyperbolic partial differential equation with variable 
coefficients, the exact type depending on <{) and <p 

X XX 

However, because flutter analysis is primarily concerned 
with the stability of small perturbations about a steady 
flow, the oscillatory component may be assumed small 
compared to the steady flow potential and therefore the 
product term be neglected, yielding. 



(l-M^)ii;^ + ipyy - M^(y+1)I(|5 




0 



(III-20) 
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IV. LINEARIZATION OF THE GOVERNING EQUATION 



The basic flutter equation, (111-20)^ is still a non- 
linear, mixed elliptic-hyperbolic partial differential 
equation with variable coefficients and difficult to solve. 
It may yet be further simplified. 

A. BASIC SOLUTION 

For M = 1, equation (III-20) becomes 




0 



(IV-1) 



Now assume 



w 



constant 



(IV-2) 



X 



0 



XX 



throughout the interblade channel . Setting 



<}) (y+1) = w(y+1) 



2 



(IV-3) 



m 



yields 



2 , 
m ip 




0 



(IV-4) 
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The solution to this equation is found in Garrick and 



Rubinow [ 5 ] 



'Mx,y) 



x-my 

~ ^ / V (s) J 

m ^ * r 



j ^ y iA(s-x) 

[0) Y (x-s) (my) '^] e ds 



X, 



for y > 0 



(IV- 5a) 



and 



x+my 



ti^(x,y) = j / v(s) k/(x-s) ^-(my) ^]e^^ 



X, 



for y < 0 



dx 



(IV-5b) 



where 



V (x) 



lira ip (x,y) 
y-0 



v(x) may be obtained directly from the tangential flow 
boundary condition, and 



0 ) 



k^(l-m^) 

4 

m 



A 



X 



Z 



— j (where this paper 

m employs the 

approximation used by 
Gorelov [4] ) 

blade leading edge. 
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Gorelov [4], has proposed a further simplification. 



Set 



z = my 



(IV-6a) 



'f(x,z) = 4^(x,y) e 



iAx 



(IV- 6b) 



Equation (IV-4) then becomes 



- '1' + w "F = 0 

XX zz 



with solution 



x-z 



'^{x,z) = - ^ / v(s) J^[ojV(x-s) ^-z^]e^^®ds 



:iV-7a) 



z > 0 



x-z 



iAs 



<'(x,z) = i / v(s) J^[ojV (x-s) ^-z^^ ® 



(IV-7b) 



z < 0 



where 



v(x) = me lim 'F (x,z) 

z-0 ^ 



v(x) is obtained from the tangential flow boundary condition. 
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For a thin body immersed in the flow, the solutions 
for y > 0 , z > 0 , and y < 0 , z < 0 apply above the body 
along left-running Mach lines, or below along right running 
Mach lines respectively. 

B. BOUNDARY CONDITIONS 

1 . Flow Tangency Condition 

The boundary condition comes from the tangential 
flow condition, (III-18) 

v(x) = [1 + ik(x - x^)] (IV-8) 

2 . Upstream Condition 

The final linearlized equation is a hyperbolic differ- 
ential equation with boundary condition 

ip(xry) = 0 (IV-9) 



when 



X - x^ < I my I 

for the solution shown in equations (IV-5) or 

'Mx,z) = 0 
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when 



X - < |z| 



for the solution shown 



in equations 



(IV-7) . 
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V. PROBLEM FORMULATION 



A. CO-ORDINATE SYSTEM 



yr 








Assume the geometry shown above. Both blades are thin 
airfoils of semichord c. All measurements are non-dimensional, 
normalized to c. The (x,y) co-ordinate system has its 
origin at the center of the reference (lower) blade. The 
(x,y^) system is centered at the middle of the adjacent 
(upper) blade. The origin of the (x,y^^) system is located 
at (0,p) in the reference system. Generalizing this conven- 
tion, the same symbols shall be used for the same quantities 
on both blades. Where discrimination is required, the quantity 
associated with the adjacent blade will be marked with 
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superscript the quantity associated with the reference 

blade will be either unsuperscripted or marked with a superscript 

Each blade is assumed to perform a small amplitude harmonic 
oscillation about its mid-chord point. Both blades are assumed to have 
identical reduced frequencies, k, and the motion. of the adjacent 
blade lags that of the reference blade by a phase angle a. 

The blades are immersed in a uniform flow from the left 
at M = 1. The objective is to determine the oscillatory 
pressure distributions and aerodynamic forces generated by 
the blades' oscillations. Cascade solidity, x = 2/p. 

B. BOUNDARY CONDITIONS 

1 . Upstream Condition 

li; = 0 whenever 

X + 1 < I my I 
and 

X + 1 < |myj^|, simultaneously 

2 . Flow Tangency Condition 
Along the reference blade 

lim tj) (x,y) = (1 + ikx) (V-2a) 

y^O ^ 

Along the adjacent blade 
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(1 + ikx)e 



(V-2b) 



lim ijj (x,Yt) 
y ^-0 1 



ia 



where 



a is the phase angle between the blades oscillations 



C. BASIC SOLUTION TECHNIQUE 

Assume that the unsteady potential, ip , may be written 
as the sum of four components 



(x,y ) 



(!)°{x,y) + + 4»^(x,yj_) + t{/^(x,y^) (V-3) 



where : 



= potential due to the reference blade 
alone, known from equation (IV-7) 

<})^ = potential due to the adjacent blade 

alone, known from equation (IV-7) 

= interference potential required to satisfy 
tangential flow condition along 
reference blade, unknown 

= interference potential required to satisfy 
tangential flow condition along 
adjacent blade, unknown. 



This total potential must satisfy the tangential flow 
condition at the plane of both the reference and adjacent 
blades. Thus 



4>y{x,y=0) + 



(x,y,=-p) + i/^°(x,y=0) + 
^1 ^ 



4'y^ (x,yi=-p) 



(V-4a) 



= (1 + ikx) 
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at the reference blade , and 



4>°(x,y=p) + <{)5 (x,y,=0) + >i^°(x,y=p) + (x,y =0) 

^ ^1 ■'* 



= (1 + ikx) e^^ 

at the adjacent blade. 

But from the unsteady potential solution for a single 
oscillating blade one has 


(V-4b) 


<Py(x,y=0) = 1 + ikx 

and 


(V-5a) 


<P^ (x,y.=0) = (1 + ikx) e^° 

Yl 1 

Thus 


(V-5b) 


<P^ (x,y,=-p) + (x,y,=-p) + >j;°(x,y=0) = 0 

^1 ^ 

along the reference blade/ and 


(V-6a) 


4>°(x/y=p) + 4'°(x,y=p) + (x,yj_=0) = 0 

along the adjacent blade. 

From equation (IV- 7) 


(V-6b) 
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(f>° (x,y) 



where 



(x,y^) 



x-my 

= ■ ^ / v° (s)JJu)l/(x-s) ^-(my) ^]e^^ 

-1 

y > 0 (V-7a) 

x+my 

= ^ / v° (s)JJojV (x-s) ^-(my) ^]e^^ 

-1 

y < 0 (V-7b) 

v°(s) = 1 + iks 

A = k/m^ 

m = (y+l)w 

k^(l-m^) 

(0 - 5 

m 

w = mean value of 4 in the channel 

x-my^^ 

= - i (s)J^[(i)V{x-s) (my^^) ^]e^^ (s-x) 

-1 

y^ > 0 (IV-8a) 

x+my^ 

= - j (s) J [ojV (x-s) (my, ) ^]e^^ 

m o 1 

-1 

y^ < 0 (IV-8b) 
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where 



(s) 



= (1 + iks)e 



10 



Henceforth attention will be restricted to the flow within 
the channel, 0£y£p , ~P£yj_f.O leaving (IV-7a) and 
(IV-8b) as the governing equations of interest. 

The two interference potentials are assumed to have 
forms identical to (IV-7a) and (IV- 8b) . 

Set 



(x,y) 



x-my 

1 r o , V -p r *\l f iX (s— x) j 

/ u (s) J [oW(x-s) -(my) le da 

m o 

-1 



y >• 0 



(V-9a) 






(x,yj_) 



x+m 

i / u^ (s) Jq [ajV (x-s) (my^)^] e^^ ^^ds 

-1 



y^ < 0 



(V-9b) 



where 



u°(s) and u^(s) are unknown functions to be 
determined so as to satisfy equations (V-6) 



Substitution of (V —7a) , (V -8b) and (V —9) into (V -6) 
yields 



O / V 

u (x) 



+ (x,y^=-p) + (x,y^=-p) 



= 0 



(V-lOa) 
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0 



(V-lOb) 



u^(x) + >|^°(x,y=p) + <|)°(x,y=p) 



Recalling Gorelov's transformation discussed in 
[4] and shown in equations (I-V-6) above, set 



Then 



, X aO, V i^x 

$ (x,z) = 4) (x,y) e 



^ 1 / N a. 1 ■ X 

4> (x,z^) = 4> 'x,yj^) e 



)iiO , V , o , , i^x 

'i' (x,z) = t(j (x,y) e 



'{'^(x,z^) = i(j^(x,y^) e^^^ 



where z = my 



Zi = my^ 



iXx 



m 



u° (x) + (x,z^=-x*) + (x,Zj^- x^) - 



= 0 



(V-lla) 



and 



where 



iXx 



u^(x) + 4>°(x,z = X*) + 'i'°(x,z=x*) = 0 (V-llb) 



m 



X* = mp 
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To employ the collocation method, assume that (x) 
and u (x) can be approximated as the sum of a set of elemen- 
tary functions f ^ so that 

n 

u°(x) :: ); 0° fj(x) (V-12a) 

j=0 

n 

u^(x) ^ f j (x) (V-12b) 

j=0 

where f .(x) = 0 when x < x*-l 
3 — * 

Note that here both u° and u^ are expressed in terms of the 
same elementary functions, f j . 

Because of the slightly supersonic nature of the problem 
observe that u*^(x) = 0 and u^(x) = 0 when x £ x* -1. 

Equations (V-12) may now be rewritten as 



x+z. 



^^^J0°f.(x) / ^0^f (s)J [wV(x-s)^-z^]e^^®ds 



x*-l 



x+z. 



= -e^^ ^ — / (1+iks) J [tJ'/ (x-s) ^-z^] e^ 

dZ, ' O i 



As 



ds 



-1 



at z^ = -x^ 



(V-13a) 
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9z 






X-2 



— f 



x*-l 



[e°f^(s)j^[wV(x-s)^ 



, iAs , 
-zje ds 



x+z 

= ^ / (l+iks)J^[ooY(x-s)^ - z^]e^^®ds 

-1 

at z = x^ (V-13b) 

where 

f .(x) = 0 for X < x+ - 1 

3 — 

This simplifies to 

x+z^ 

eiAxrgOf (x) ^ yQ^{_A_ I f.(s)J [wV(x-s ) ^-z? ] e^^^ds } 
'id ^ D jo 1 



x+z^ 

/ (1+iks) [(jS^(x-s) ^-z^] e^^^ds 



at z^ = -X* (V-14a) 

and 

x-z 

e^^^]»9^f j (x)-^0°{-^ / f j (s) J^[u)V (x-s) ^-z^]e^^®ds 

x*-l 

X-Z 

= 9 / (l+iks)j [ 0 )V(x-s)^-z^]e^^®ds (V-14b) 

d Z O 

-1 
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at z 



= X* where 



f ..(x) = 0 for X < x*-l. 
1 — * 



Performing the indicated differentiation yields 



i^x 






x-x. 



3 3 



fAs) 



x*-l 



[(ijY (x-s) ^-x^] (jJx*e^^^ds 

T 2 2 

V(x-s) -X* 



+ f .Cx-x.)e^^<='-=‘*> ) 



x-x. 



la 



r r . "★'"1 

e {- / (1+iks) 

-1 



cox.Jt [ojV (x-s) ^-x?] e^^^ds 



-/ (x-s) ^-x* 



- [1+ik (x-x*) ] e^^ } 



(V-15a) 



and 



x-x^ 



t le°{ / £j(s) 



J. [wV(x-s)^-x*l(^x*e' 'ds 



iAs 



x*-l 



V(x-s) 



+ f . (x-x*) e 



iA (x-x*) j 



x-x. 



/ (1+iks)' 



oox 



-1 



(x-s) ^-x*e^ ^ds 
V (x-s) ^-x* 



- [1+ik (x-x^) ] e 

(V-15b) 



iA (x-x. 
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Gorelov's formulation, equations [2.6, 2.7, 2.8, and 2.9] 
of [4] , can be obtained directly from equations (V-15) by 
substituting 



fj (x) = cos jo - cos jo* 

0 r' O • 

V = ) U- cos 30 

j=0 J 

1 ^1 

V = I cos jn 

j=0 ^ 

where : 

0 = cos (-x) 

n., = cos (1-x*) 

In comparing the two systems care must be taken to note the 
differing symbols and coordinate systems. The corresponding 
quantities are: 



Here 


in [4] 


e°, ej 


V , V. 

oa la 


O 1 

p . , y . 
J J 


0 , 9 ■ 

oa la 


0 = cos ^ (-x) 


n = cos ^(1-x) 


^ = cos ^(1-x*) 


0 * = cos ^(1-x* 


z , z^ 


1 — 1 
>1 

>1 


j 


a 


a 


4' 
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Here -1 £ x £ 1 ; in [4] 0 _< x _< 2 . This transformation 
accounts for the differing definitions of n* Making the 
substitutions results in the system 



n 

I { 0 °[cos jn - cos jn*] 

j =0 

(l-6io)cos jn*le^^®ds 



x-x 



+ 0 - 







[(J 0 ^l(x-s) [cos jn - 



^Ir ■— /n X iA(X-X.)x 

+ 9t[cos* jn - cos jn* ] e * } 



n 



x-x. 



/ 

j=0 -1 






o 



V(x-s) 



2, i^s, 

z^] cos Jne ds 



1 iA (x-x*) 1 

IJj cos jne * 



(V-16a) 



and 



at 23 ^ = -X* 
X > x*-l 
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^{9^ [cos jn - (1-6 .) cos jn+] 
J ^ 



x-x , 



+ 9° 






/ ^ Jq [oj'V (x-s) ^-z^] [cos jn~ (1-6^^ 



1 iAs, 

cosjn^Je ds 



+ 9°[cos jn. - (l-6^^)cos jn*J } 

n x-x* 

I Jq [oiV (x-s) ^-z^] cos jne^^^ds 

j=0 -1 



where : 



- u° cos 



(V-16b) 



n = arc cos (-x) 

n* = arc cos (1-x*) 

n = arc cos (-s) 

✓V 

n = arc cos (-x+x*) 

6 - . = Dirac 6 function 

<?3 



1 for j = 0 

0 for 3^0 



Given the change in coordinates and notation, this system 
is equivalent to that shown in [4] . 
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This was the system programmed for computer solution. 



Because the function, fj(x), is unaffected by the differen- 
tiation with respect to y (or z) the exact form used need not be 
specified, so that the system shown in (V-15) may be programmed 
with f undetermined. A subroutine may be written to return 
the function desired and the remaining program left perfectly 
general. In the program developed with this thesis both the 
Gorelov functions shown above and the Legendre polynomials 
were employed. All the integrals may now be evaluated at 
n+1 points, on both blades in mp-1 < x^ < 1 and the 

resulting linear system solved for 9° and 0^, j=l,2, . . . ,n+l. 

The interference potentials may be constructed by taking 



X- z 





1 



z > 0 



(V-17a) 



x+z 



1 




m 



1 



/ 




-1 




(V-17b) 



where 



f j(x) = 0 , for all 



X < x*-l 



Once the potentials have been calculated as outlined 
above, the surface pressure may be calculated using the 
relationship 
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(V-18a) 



Cp = + ikil^) 

= -2[^^ + i(k-A)H'] (V-18b) 

Because all the plates are assumed to be in steady oscillation 
with uniform phase shift, a, between neighboring plates, 
then 

v^(x) = v°(x)e^^ , u^(x) = u°(x)e^'^ , ''F (x,y) = -’F{x,-y) 

in this case 



1 

“"LA 

Cl = 2 / [H- (x,+0) + i(k-A)'F(x,+0) ]e ‘^'^^dx (V-19) 

Ct X 

-1 



where 



'FUj_.+0) 



X 

^ / [(l+iks)+u°{s)]J^[uj(x-s)]e^^^^^ds 

-1 






^ . - __ ^ 

^ / [(l+iks)+u°(s)]J [o)V(x-s)^-(x*)^]e^^^®^ds 

m 

-1 



where x* = mp. 



(V-20) 



Results from this approach, in the form of values of 
for k = 0.1, at various values of w are presented in the results for 
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approximations based both on Gorelov's formulation, and on 
the Legendre polynomials . 



D. COLLOCATION SOLUTION OF THE POTENTIAL EQUATION EXPANDED 

FOR SMALL k, 

In order to provide a partially independent check of the 
results of the main program, the Gorelov function represen- 
tation of the collocation solution was expanded for small 
k, and solved at two collocation points, n = 2. The resulting 
potentials, and partial derivatives with respect to x and 
y were then used to replace the corresponding niomerical 
routines in the main program. The output resulting from the 
approximations were compared with the purely numerical results 
obtained from the computer program. 

1. Solution For The Unknown Potential Coefficients 



The basic system of linear equations used to determine 



the unknown coefficients is 



m 



1 iAx o 
— e u 




1 



+ 

z 



1 



0 , z = 0, ‘^1 ~ ~ 



(V-21a) 



1 iAx 1 
— e u 



+ = 0 
z z 




(V-21b) 



m 



where 



0 = u°(x) = u^ (x) when x £ -1+x* 
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otherwise 



n 

= I 9^ (cos jn - cos jn*.) + 9° 

j=l 

n 

1 r 1 1 

u = ^ 9j (cos jn - cos jn*) + 9^ 

j=l 

where 

n = arc cos (-x) 
n* = arc cos (1- x*) . 

Thus, for n = 2, the system becomes 

eiAX{gO 9° (cos n ~ cos n*)l 

^ X+2, , , J ^ _ 

0 

(V-22a) 



3z. 



/ ^ [o)V(x-s)^ - z^je^'^^ds 



-1 



x+z • 



3z. 



/ (s) [a)V(x-s) ^-z^] e^^^ds = 



x*-l 



Zl = -X* 



along the reference blade and 
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x-z 



eiAX{gl 0^(cos n - cos n*) } 
O i ^ 



/ v° (s) [o)V(x-s) ^-z^] 

-1 



x-z 



3z 



/ u°(s)J [col/cx-s)^ - Z^le^^^ds . 



x*-l 



Z = Xj 



along the adjacent blade 
o 



where 



u'^(s) = + 6° (cos n - cos n*) 



(V-22b) 



u^(s) = 9^ + 9^ (cos n - cos n*) 

o i 



V (s) = 1 + iks 



(s) = (1 + iks)e^‘^ 



For k sufficiently small, this system may be further 
simplified by the following approximations 



J lo) (x-s) ^-z^] = 1 - O(CJ^) 

o 



(y-23a) 



J [coV (x-s) ^-z^] = 1 - 0 (co ) ~ 1 



(V-23b) 



e^^^ :: 1 + iAX - Q(A^X^) :: 1 + iXX 



^ 1 + iAs - 0(A^S^)_ :: 1 + iAS 



2 2 
where 0 (co ) means "of the order of magnitude of cj " 



-l£x£l, x-x* 

^ /_2 2 _ k^(l-m^) 

A = K/m , 0 ) = J 



m 
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The interference source distributions may be replaced by 



u°(s) = 9 q + 0°(-s + x*-l) 

u^(s) = 0^ + 0^(-s + x*-l) . 

O 1 * 

If higher order terms are neglected, the result is a system 
linear in k and A 

x+z, 

3 ia 

(1+iAx) [0° + 0°(-x-l+x*)] +11 — / (1+iks) (l+iAs)ds 

^ -1 



+ 




x+z^ 

/ [0^ + 0^^ (-S-1+X*) ] (1+iAs) ds 

x*-l 



z 



1 



-X 



★ 



0 

(V-24a) 



x-z 



(l+iXx) [0Q + 0]^(-x-l+x*) ] 



- — / (1-iks) (l+iAs)ds 

o Z 

-1 



x-z 

3z ^ o 

X*-l 



0°(-s-l+x*) ] (l+iXs)ds 



= 0 



Z = X. 



(V-24b) 



Product terms containing (kX) = ^ may be neglected as 

m 

of higher order in k, yielding 
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x+z. 



la 



(1+iAx) [9^ + 9j^ (-X-1+X*) ] + — I [1+i (A+k) s] ds 



3z. 



-1 



x+z . 



/ nej + e];(-s-i+x 



.)] 






+ iAs[9^ + 9^ (-s-l+x^ ) ] }ds = 0 

= -X* (V-25a) 

x-z 

(l+iAx)[9^ + 9^(-x-l+x*)] ~ / [1+i (A+k) s]ds 

-1 

2 

{[9° + 9°(-s-l+x*)] + iAs[9° + 9° (-s-l+x^ ) ] }ds 

z = X* (V-25b) 

Evaluating the indicated derivatives yields 

(1+iAx) [9° + 9?(-x-l+x*)] + [ 1+i (k+A ) (x-x* ) ] e^*^ 
o 1 * ^ 

+ [9^ + 9^- (2x*-x-l) ]+iA (x-x*) [9^+9?- (2x*-x-l) ] } = 0 
O 1 ^ O 1 ^ 

(V-26a) 
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(1+iAx) [0^+0 (-X-1+X*) ] + [1+i (A+k) (x-x*) ] 

+ { [0°+0° (2x*-x-l) J + iA (x-x*) [0°+e°(.2x*-x-l) ] } = 



Thus : 

0°(l+iAx) + 0°[(-x-l+x*) + iAx(-x-l+x*) ] 

+ 0^[1 iA(x-x*)] + 0]^[ (2x*-x-l) + iA(x-x*)(2x* 

= + i(k+A) (x-x*) J 

0q[ 1 + iAx] + 0^[(-x+l-x*) + iAx (-X+1-X*) ] 

+ 0° [1+iAx (x-x*) ] + 0° [ (2x*-x-l) +iA (x-x*) (2x*-: 
= -[1 + i(A+k) (x-x*) ] 

This system may be solved at two points, x^, and X 2 / for 
0°, 0°, 0^, and d^. 



0 

(V-26b) 

-x-D] 

(V-27a) 

- 1 )] 

(V-27b) 
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2. Calculation Of The Potential 



The potential is given by 



^ (x,y) 



'i'(x,z)e"^^^ 



(V-28) 



where 



X 

'i'(x,z) = - ^ / [v° (s)+u°(s) ] J^[toV(x-s) ^]e^^®ds 

-1 



x-x* 

+ ^ / [v^(s)+u^(s) ] J^[coV(x-s)^-x*Je^^®ds. 

-1 



u°(s) = 



(s) 



0 for all s £ x^-1 . 



Thus 



X 

i'(x,z) = “ ^ I v° (s) lo)V(x-s) e^^^ds 

-1 



X 

- i / u°(s)J^[u)V(x-s)^]e^^®ds 
x*-l 

x-x* 

+ ± f v^(s)J^[ojV(x-s)^-x*]e^^^ds 

-1 



x-x. 



4 / 



u^ (s) J, 



[o)V(x-s) ^-x?^ ® 



iXs 



ds (V-29) 



x*-l 
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Making the same small frequency approximations as 



in the previous section yields 



X 

'{'(x,z) = - i j v°(s) (l+iAs)ds 

/ lU 

-1 



1 

m 



x 

/ u° (s) (1+iAs) ds 

x*-l 




x-x^ 



(s) 



(1+iAs) ds 



-1 



+ — / u^ (s) (l+iAs)ds (V-30) 

x*-l 



From the general formulation 



' 1 ' = 



(V-31) 



Thus , along the reference blade 



x 

-m'J>*^(x,z=0) = / v°(s) ( 1+iAs )ds 

-1 



X 

/ (1+iks) (l+iAs)ds 

-1 



X 

/ [1+i (k+A) s]ds 

-1 



2 X 

[s+i(k+A) VJ 

-1 



(V-32) 



X + 



. (k+A) 

1-^ — + 



1 - i( 



k+A ' 
2 ’ 
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(V-33) 



$°(x,z=0) = - i[(i+x) + i(iS±i) (x^-i)] 



By inspection 



$ (x,z^=-x*) 



e 



.10 



^^{(1+x-x^) + i(lS±i) [ (x-x*) ^-1] } (V-34) 



X 

(x, z=0) = / u° (s) (1+iAs) ds 

x*-l 



(V-35) 



X 

/ [0° + 6° (-s+x^-1) ] (1+iAs) ds 

x*-l 



X 

/ 6°(l+iAs) 

x*-l 



+ 0°(-s+x*-l) (l+iAs)ds 



= 0°((x+l-x*) 

o ^ 



-^(x^-x*+2x*-l) ] 



X 

+ / 0° [-x+x*-l+iA (-s^+sx*-s) ] ds: 

x*-l 



= 0 [(x+l-x*) 

o * 



iA . 2 2 , « 1 \ 1 

-^(x -x*+2x*-l) ] 



+ 0 



0 

1 



{ (-^-^-sx^-s) +iA (- 



3 2 

s , S x^ 
" 2 




X 

} 

x*-l 

(V-36) 
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9°t (x+l-x*) +^(x^-x^+2x*-l) ] 



o. x2 (x*-l)" 

+ 9^{ (-Y“i-x*-x) - [ ^ t-x* (x*-l) - (x*-l) ] 



3 X X* ^2_ (x*-ll (x*-l)2 (x*-l)2 



, . , r / X X , , 

+ ^ ^)4 



3 3 - ^]} 



9° [ (x+l-x*) +^(x^-x*+2x*-l) ] 



+ 9° -^x(x^-l) + ^^-^~^^ (x*-l) 



+ iA [- -^ + X — 2 



2 (x^.-l) , ^ _ 



X* + X* - y 



ill + x^ - ^ + — 
2 2 2 



- X* + ^] } 



9°[(x+l-x*) + ^(x^ 
o ^ z 



X* + 2X* 



- 1 )] 



+ 9°{- ^ + x(x*-l) 






. , , x"^ x"^(x*-l) _ ^ ^ _ -E* -I- — ] } 

+ iA[- ^ r 2 2 6^ ^ 



60 



~ m^®o ^ X* + 2x*-l)] 



+ 0°{- ^ + x(x*-l) - 



3 X (x*-l) (jj _i)3 

. iA[-^ + — ^ ^L]}} 



(V-37) 



mt{) (x,Zj^=-x^) may be evaluted by substituting x-x* for x in 
the expression for m'i'° and exchanging 9^ and 9^ for 9° and 
9? 



ra’i' (x,Zt=-x^) = 



9^{ [ (x-x*) +1-X*] + ^[(x-x*)^ ~ X* + 2x* - 1] } 



t 9l(-Ul|ai + (x-x.) (X.-1) - 



. ;w (X-X*)^ , X^(x*-1) (X*-1)^,1 

+ lAl- 2 + 2 6 



(V-38) 



9^ (x+l-2x*) + (x^- 2 xx*+X 2 ) - x* + 2x* - 1] } 



+ el(- (x^-2xx,-i-x ll ^ ^ 1) 



x~ 2. 
3 



2, (x*-l) 



+ iA[- ^ + x"x* - XX* + ^ + (x^-2xx*+x*) 



(x^-1) 



-]} 
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1 ] 



= 9j{(x+l-2x*) + - 2xx* + 2x* - 



X ^ ^(^-1) + |(l-3x*) (x*-l) 



+ iA [- ^ + + x(- 



+ ^ + |]} 



4» (x,z,=-x*) 



Irnl 



iA r 2 



= -{e.{ (X+1-2X*) + V[x - 2xx* + 2x* - 1] 



+ e]^{[- ^ + x(2x*-l) - XX* + 2x* - |] 



+ iA (- ^ + + xi::^4^ + 5^ 

3 2 2 6 



i])} 



(V-39) 



A comparison of the results for the 


full program and 


the 


approximation is given below for k ^ 


= 0.01, w = 0.05, 


a = p. 


2 - 3 

n = 2, yielding oj - 6.1x10 , ^ ~ 


0.083 




Full program 


Approx 




X = .1285 


4) -5.363, .327i 


-5.379, 


.5491 


4) -8.634 , . 361i 


-8.66, 


.3791 


X 

X = .5643 


4) -9.754, .8091 


-9.804, 


l.OOOi 


d) -12.234, .692i 


-12.272, 


.76081 


' X 


C. +31.658,-3.10321 


31.7019, 


-3.21651 


A/ 

a 


a 





62 



VI . RESULTS 



The collocation method was used to solve the partial 
<^if f s^sntial equation resulting from the Gorelov approximation 
of transonic potential flow in an unstaggered cascade. The system 
was solved using both the spanning functions proposed by Gorelov in 
[4] , resulting in the equations (V— 16 ); and the. Legendre poly- 
nomials, resulting in equations (V-15) with f^ replaced by 
the Legendre polynomial, P.. The resulting values of C a 
for k=.l, X =1, a=l and seven collocation points on 
each blade are presented in figures VI-2, VI-3, and VI-4. 

Figure VI -1 presents a diagram which is useful in 
commenting on the other results. This shows the location of 
the collocation points and first three interference reflec- 
tions as a function of w expressed as a percentage of that 
portion of the chord subject to reflection. The collocation 
points are equally spaced throughout this interval, 12.5% 
from the leading edge of the interference zone, 12.5% 
between each pair of points and 12.5% from the blade trailing 
edge. The independent variable, w, is plotted vertically 
so that the dependent variable, percent of chord subject 
to interference, may be more conveniently visualized along 
the blade. (The curves are not precisely linear, but are 
very nearly so in the range shown.) 

Figure VI -2 shows the calculated with k = 0.0 in 

comparison with the results obtained from Ackeret theory. 
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Agreement is good where there is no reflection and the 
portion of the blade subject to interference is affected 
by a constant interference potential, w ^ 0.11. Throughout 
the rest of the curve the results calculated here oscillate 
above and below the theoretical values. This appears to be 
due to the discrete nature of the approximation used in the 
collocation method. Rarely is the fraction of the chord 
subject to interference reflection equal to the fraction of 
the collocation points which feel it. Where the collocation 
point fraction lags, as near w = 0.6, the collocation results 
are lower than those due to Ackeret theory. When the 
collocation point fraction leads, as it does for w _< .04 
and briefly for w = .08, the the collocation results are 
higher than those due to Ackeret theory. The fault appears 
to be an intrinsic feature of the small number of points 
sampled. This results in a set of coefficients similar to 
those which would be obtained from a generalized Fourier 
series based on the integration of the Taylor series 
expansion about each point. This obviously cannot be a good 
approximation when both the function and its derivative are 
discontinuous at the reflections. 

Figures VI -3 and VI -4 show the results of using 
Legendre polynomials and Gorelov's functions as spanning 
functions. The results for both formulations are identical. 
Gorelov's results are presented for comparison. Agreement 
is good for w > 0.05 except for an anomalous point, marked A 
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It is believed that this anomaly is due to the location of 
the first reflection just ahead of the last collocation 
point (cf . "A" on Figure VI -1) . This will yield a very 
small contribution from the reflection potential to the 
linear system from which the collocation points are 
determined. The resulting system will have a large dynamic 
range and may be ill-conditioned. 

The discrepancy between these results , and those in 

[3] for w < 0.5 is still unexplained, as is the outlying 
value for w = 0.5. 

The discontinuities in the imaginary results are believed 
to be due primarily to the reflection/collocation inter- 
action explained above. 
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Percent of chord subject to interference 

FIGURE VI-1. Location of Reflections and Collocation 
Points Shown as Percent of Chord Subject 
to Interference 
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FIGURE VI-2. Comparison of CHq^-vs-’w to that Obtained 
from Ackeret Theory 
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FIGURE VI-3. Plot of Cjla'VS-w, Legendre Polynomials 
k = 0.1, X = 1.0, a = Ti, n = 7 
compared with Gorelov's results 
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VII . RECOMMENDATIONS 



There are two recommendations to be made about the 
techniques used in the collocation method, and a new area 
in which it might be employed. 

The program developed in the course of writing this 
thesis employs adaptive Simpson's integration to calculate 
the elements of a completely determined system. This system 
is solved to provide the coefficients of the spanning functions. 
Two improvements may be made; 

1. The Simpson's integration scheme may be replaced by 

a Gaussian integrator. Experience has shown that several 
thousand function evaluations are required by the Simpson's 
integration routine when is to be evaluated for small 

w. This entails large amounts of computer time and leads 

to increased accumulations of numerical error. Use of Gaussian 
integration would probably improve both of these characteristics 
with little loss of accuracy. 

2. The present program treats a completely determined sys- 
tem of dimension 2n+l by 2n+l, and then solves that system to 
find the collocation coefficients. This procedure has worked 
satisfactorily in this thesis, but may not work as well at 
higher frequencies where the final linear system of equations 
may be ill-conditioned. As an alternative, it is recommended 
that the boundary conditions be applied at more than n points, 
say m points, where m is twice or three times as many points, 
and that the least squares technique be used to determine the 
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the collocation coefficients which give the minimum square 
error over-all. This may be thought of as "sampling more 
data" in order to get more information about the unknown 
function. The present program could be easily modified in 
this regard by replacing the spanning function matrix, QlZINT, 
by a new matrix of the form 



QlZINT' = x'^X 



where X is the new m by n+1 (m > n+1) matrix, and replacing 
the present right-hand-side vector, QICOF with 



QICOF ' = x'^Y 



where Y is the new m by 1 right-hand-side vector ., An 
alternative would be to employ a prepackaged statistical 
linear regression routine after either modifying the routine 
to accept complex data, or transforming the present system 
into a larger system of real numbers only. 

The new area in which the collocation method might be 
employed is the calculation of the potential flow about a 
staggered cascade. The method coud be employed to calculate 
both the potential in the channel and above the upper blade. 
The program presented has been designed to enable the 
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calculation of flow within the channel of a staggered 
cascade. Unfortunately, there was not enough time to extend 
the study to this case. 
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APPENDIX A 



PROGRAM DESCRIPTION 

This section describes the computer program used to 
calculate the interference solution to the Gorelov lineariza~ 
tion for unsteady transonic flow in a channel. The program 
written in IBM Fortran IV with the basic structure outlined 
by Stevens [5] . The basic points are; 

1. Organization of the program into small subroutines, 
each of which performs a specific task. 

2. Transmission data to and from subroutines via a 
formal parameter argument list. No common statements 
are used. 

The end objective is code which is both easy to modify and 
maintain . 

Each subroutine is designed with optional diagnostic 
printing of its input and output. This is controlled by the 
parameter IPT. The diagnostic output is printed (only) if 
IPT > 0. Each routine accepts IPT, sets lOT = IPT - 1, and 
then passes lOT as the print parameter to routines it calls. 
By this method, diagnostic output can be "cascaded" to any 
desired level. Large initial values of IPT should be avoided 
because of the spectacular amount of output which can be 
generated by the double integrals within QIDCOF. 

1. Main Program ; including subroutines READ and ABSA. 

The basic structure is given in Figure A-1. MAIN 
calls READ to read input data and then ABSA to calculate the 
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MAIN 



READ 



ABSA 



Q$ONE 



QIZINT 

— QICABP 
— QICRBP 
— QICAIP 
— QlCRIP 



QICOEF QIPPOT QIDCOF 



LEQTIC 




QIDABP 

QIDRBP 

QIDRIP 

QIDAIP 



QIDABP 

QIPRBD 

QIPAID 

QIPRIP 



Figure A-1. Program Hierarchy 



collocation points. The version of ABSA shown evenly spaces 
the collocation points across that portion of the blade 
siabject to interference. ABSA may be easily replaced if 
different point spacing is desired, or if additional points 
are to be added for an overdetermined system and least 
squares approximation. 

2. Q$ONE This sxibroutine controls the actual potential 
calculation. It performs no calculation itself, but calls 
subordinate subroutines where the calculations are actually 
performed. The calling hierarchy is shown in Figure A-1. 

3. QIZINT This siibroutine calculates the linear equa- 
tion system arising from the boundary conditions. Hierarchy 
is shown in Figure A- 2. 

The matrix output is carried through QIINT. QIZINT 
calls the following subprograms 

a. QlCRBP returns the value of 

b . QICABP returns the value of (J) ^ 

^1 

c. QICRIP returns the values of \l>° 

X-2 

^ (x-s) ^-z^Jds 

dz ■'Jo 

-1+x* 

where f^ is one of the set of i elementary 
functions, j=l,n 

1 

d. QICAIP returns the values of ip 
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Figure 



x+z 




3z 



/ 



1 



-1+OFFSET+x* 



where f j is one of the set of elementary 

functions. OFFSET is a parameter 
included to facilitate program 
conversion to a staggered cascade 

e. PLGNDR is the subprogram which returns f . (x) , 



the elementary spanning function. No other routine contains 
explicit reference to the spanning function. This facili- 
tates easy replacement of the spanning functions should this 
be desired. 



QICWXM computes end-points and then calls QlCZXM, a complex 
integration routine based on SIMP by Shampine and Allen [6] . 
QlCAIP and QlCRIP call QlCWSM and QlCZSM to- perform the 
integration. QIZINT passes the constant matrix to QICOEF 
in the array QIINT and the right-hand-side vector in the 
array QICOF. 

4. QICOEF This subroutine employs the IMSL routine 
LEQTIC to solve the linear system received from QIZINT. The 
resulting coefficients are QIABCF for the adjacent blade and 
QIRBCF for the reference blade. LEQT2C, the high precision 
complex IMSL routine may be directly substituted for LEQTIC. 
QICOEF may be rewritten to employ the generalized inverse 
required for least squares approximation 



QICRBP and QICABP in turn call QICWXM and QlCZXM. 
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A = (x x) X y where x = QIINT 

y = QICOF 

after first performing the multiplication necessary to replace 
QIINT and QICOF with the proper matrix products in the call 
to LEQTIC. 

5. QIPPOT This subroutine calculates the potential, 

(J), and <j) , at each collocation point along the reference 
blade but only if QIPPOT receives a value of IPT > 0, re- 
quiring IPT ^ 2 on input to the main program. If IPT £ 0, 
then the subroutine is exited before any calculations are 
performed. This subroutine is most useful for debugging 
QIZINT and QICOEF. QIPPOT calls QIPABP, QIPRBP, QIPRIP, 
QIPAIP,. QIDABP, QIDRBP, QIDRIP , and QlDAIP , all of which will 
be described in the next section. 

6. QIDCOF This subroutine calculates the dimension- 
less coefficients of lift and moment; a' *“^a‘ 

internal hierarchy is shown in Figure A-3. 

a. QIDCL calculates the nondimens ional complex 
coefficient Co 

b. QIDCM calculates the nondimens ional complex 
coefficient 

c. QIPRBP and QIPABP calculate the potentials due 

to the reference and adjacent blades respectively. 
QIPWXM and QlPZXM are called to perform the 
actual integration. 
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d. QIPRIP and QlPAIP calculate the interference 
potentials along the reference and adjacent 
blades. QIPZSM is called to perform the 
integration. 

e. QIDRBP, QlDABP, QlDRIP, and QlDAIP correspond 
exactly to subroutines above except that the 
value returned is the partial derivative of the 
potential with respect to X. QIDWXM, QIDZSM, 
and QlDXSM perform the co-reponding integrals . 

6. Program Listing The program listing shown below 
incorporates the Legendre functions as spanning functions . 
Listings for a subroutine employing Gorelov's spanning 
function and the linear approximation program follow. 
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7 . Gorelov Spanning Function . 

The subprogram for function used by Gorelov is 
shown below. 



RE^L FU^CTICN PLGNOR+8 (X, OR tN ) 
INFLICIT RE4L*8( A-I-,C-Z) 
IF(N.SQ.O) GOTO ICO 
ET/ = CARCCS(-X) 

ElASTR = DARCOSd.OOO - DR) 

FN = CFLCAT(N) 



IOC 



FLG^DR 

RETURN 

PLGNCR 

RETURN 

END 



= DCCS(FN*ETA)-OCQS( FN*ETASTR» 

1 .000 
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8 . Linear Expansion Program 

The program based on the linear expansion for 
small k is shown below. 
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